Detection of supernova neutrinos on the Earth for large #13 
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Supernova (SN) neutrinos detected on the Earth are subject to the shock wave effects, the 
C**) ■ Mikheyev-Smirnov-Wolfenstein (MSW) effects, the neutrino collective effects and the Earth mat- 

■ ter effects. Considering the recent experimental result about the large mixing angle 613 (^ 8.8°) 
provided by the Daya Bay Collaboration and applying the available knowledge for the neutrino 

■ conversion probability in the high resonance of SN, Ph, which is in the form of hypergeometric 
function in the case of large #13, we deduce the expression of Ph taking into account the shock wave 
effects. It is found that Ph is not zero in a certain range of time due to the shock wave effects. After 
considering all the four physical effects and scanning relevant parameters, we calculate the event 

i-Cj , numbers of SN neutrinos detected at the Daya Bay experiment. From the numerical results, it is 

found that the behaviors of neutrino event numbers detected on the Earth depend on the neutrino 
, mass hierarchy and neutrino spectrum parameters including the temperature T a , the dimensionless 

pinching parameter 77c or fi a (where a refers to neutrino flavor), the average energy (E a ), and the 
SN neutrino luminosities L a . We also compare the results of two parametrization methods for the 
neutrino energy distributions and give the ranges of SN neutrino event numbers that will be detected 
at the Daya Bay experiment. 

O . PACS numbers: 14.60. Pq, 13.15. +g, 25.30.Pt, 26.30.-k 
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In the universe, core-collapse supernovas (SNs) are among the most energetic explosions [H-Q- They not only 
mark the catastrophic end of some stars, which turn into neutron stars or black holes after explosions, but also are 
responsible for the richness of heavy elements Q-Q- SN1987A has attracted worldwide interests and has been studied 
extensively since it came into our sight several decades ago Qfl- During the explosions of the type II SN, most of 
the binding energy is released as neutrinos, which are very useful for acquiring information about intrinsic properties 
and the explosion dynamics of SN 

For the past few decades, in most theoretical models, it has been believed that the neutrino mixing angle #13 is 
smaller than 3°, or even smaller than 1.5°. In Refs. [9(|l0j], the authors studied possible methods to measure this 



neutrino mixing angle while #13 < 3°. However, in recent years, some new experimental results indicated a large #13 
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This year the Daya Bay experiment measured the value of 6*13 to 5.2a 
±0.. 
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171 ] . which is much larger than 3°. Table I shows a summary 



(by large 13 we mean 13 ^ 9°) [ll[-|l 
accuracy and obtained the result #13 
of the recent experimental results about #13. In this paper, we will study detection of SN neutrinos on the Earth in 
the case of large #13. 



TABLE I: Summary of experimental results about 613. Numbers with (without) brackets are for normal (inverted) mass 
hierarchy. 



Experiment 


Result 


Accuracy(Sensitivity) 


Year 


Reference 


Double Chooz 


0.03 < sin 2 26»i3 < 0.19 


0.03 


2006 


M 


MINOS 


< sin 2 20i3 < 0.12(0.19) 


la 


2009 


[12] 


T2K 


0.03(0.04) < sin 2 26»i3 < 0.28(0.34) 


2.5cr 


2010 


[13] 


MINOS 


2 sin 2 26*23 sin 2 20 13 < 0.12(0.20) 


90% 


2011 


[14] 


Daya Bay 


sin 2 26»i3 = 0.092 ± 0.016(stai) ± O.OOS(sysi) 


5.2cr 


2012 


[15] 


RENO 


sin 2 26*13 = 0.086 ± 0.041 (siot) ± 0.03(sysi) 


4.9cr 


2012 


[16] 



SN neutrinos arc produced from the core collapse of SN and propagate outward to the surface of SN. Then they 
travel a long cosmic distance to reach the detector on the Earth. During this process, they pass through the SN 
matter and the Earth matter. While they pro pagate , SN neutrinos are subject to the Mikheyev-Smirnov-Wolfenstein 
Q Q , the shock wave effects 



2l| , the neutrino collective effects 



22 1 , and the Earth matter effects 



(MSW]effects 

[9 1 10 1 23j . Different from the case in vacuum, the behavior of neutrino oscillation changes while neutrinos propagate in 



matter. The neutrino matter effects, due to the interaction between matter and neutrinos, was found by Wolfenstein, 
Mikhcyev and Smirnov, and was named as MSW effects. Inside the SN, the large mixing angle solution of neutrinos 
results in the crossing probability at the low resonance region Pj, ~ 0. Therefore, we only need to consider the 



crossing possibility at the high resonance region, P, 



'h hi 



In addition, the shock wave effects may change the 



density distribution of SN and the position of the high resonance. The expression of Ph for large #13 needs to be 



developed 
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be considered 



26] 



19 251 . When studying the Earth matter effects, the realistic density distribution of the Earth needs to 



27j . By scanning the ranges of relevant parameters appearing in the neutrino spectra, it is possible 



to obtain the maximum and minimum event numbers of SN neutrinos detected on the Earth. There are two different 
parametrization forms for the neutrino energy distribution [l| j^. We will compare the results from these two models. 

This paper is organized as follows. In Section II, we present a very brief overview of SN explosions and the production 
of SN neutrinos. In Section III, the four physical effects on the detection of SN neutrinos including the MSW effects, 
the shock wave effects, the neutrino collective effects, and the Earth matter effects are discussed respectively. In this 
section, the expression of Ph in the high resonance region for large 6*13 is obtained. In Section IV, using the latest 
result of the Daya Bay experiment, 6*13 = 8.8° ± 0.8°, we take into account the four physical effects and calculate the 
event numbers of SN neutrinos detected on the Earth. The results of two different parametrization methods for the 
neutrino energy distribution arc compared and the effects of relevant SN neutrino parameters arc discussed in detail. 
Finally, in Section V we give the summary and discussions. 
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II. SN EXPLOSIONS AND SN NEUTRINO SPECTRA 



According to the presence or absence of hydrogen lines in their spectra, supernovas can be classified into two 
types, type I and type II. The so-called type II supernovas have hydrogen lines in their spectra, and are core-collapse 
supernovas. In this paper, we only pay attention to the type II supernovas, which are one main source of neutrinos 
in the universe j^ . The explosion process of core-collapse supernova can be divided into several phases, and more 
details about the scenario of explosion can be found in Ref. [l| . 

An SN explosion approximately releases a total energy of Eb — 3 x 10 53 erg, about 99% of which is radiated away 
as SN neutrinos 



given by 



m 



301 ] . The relation between the total SN energy and the luminosity of different flavor neutrinos is 



£«,.(*) + £».(*)+£*.(*) = — e 

T 



-t/r 



(1) 



where v x represent u^, u T , and v r . The luminosity flux of the SN neutrinos L a (a = v e ,i> e ,i/ x ) decays in time as 
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32] 



33 



L a (t) = L° a e-^. 



(2) 



The range of r was obtained by fitting the experimental data of SN1987A: r = 1.74 — 4.19s [30I] [32I] . 

In general, the SN neutrino spectra are parameterized in two forms to match the result of Monte Carlo simulations. 
One is known as the Fermi-Dirac distribution 34j . 



F(°)(E) 



E 2 



F„ T 4 e ( E / T °-')° 



where F Va is defined as: 



(3) 



-dx, 



exp (x - r) a ) + 1 

and E, r) a , T a are the neutrino energy, the pinching parameter of the spectrum, and the neutrino temperature, respec- 



tively. The spectra obtained from numerical simulations can be well fitted by choosing 
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T Vc = 3 - AMeV, T Vc = 5 - 6M eV, T„ a = 7- 9MeV, 
Vv* » 3 - 5, ^ M 2.0 - 2.5, t?„. « - 2. 



(4) 



In addition, the luminosity ratios of different flavor neutrinos, which will be shown later to play an important role in 
the calculation of event numbers, are as follows: 



|^ = 0.5-2, ^=0.5-2. 



(5) 
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The other parametrization form of SN neutrino fluxes given by the Garching group can be expressed as 

where (E a ) is the average energy of neutrino and j3 a is the dimensionlcss pinching parameter. For different neutrinos, 
their ranges of values are obtained as 



(E Ve ) = 12- 15MeV, (E Pe ) = 12 - IhMeV, (E Vx ) = 15 - 18MeV, (7) 

(3 a =3.5-6. 

The ranges of luminosity ratios for different flavor neutrinos arc the following in this model: 

-^=0.5-0.8, -^=0.5-0.8. 

(8) 

III. PHYSICAL EFFECTS ON THE DETECTION OF SN NEUTRINOS 

When neutrinos propagate outward to the surface of SN, they can be subject to the SN shock wave effects, the 
MSW effects, and the collective effects. Before arriving at the detectors, they travel through the Earth matter and 
are affected by the Earth matter effects. In this section, we shall consider all the above four physical effects on SN 
neutrinos. 



A. MSW effects and conversion probability 

The MSW effects are caused by neutrino interactions with matter, which are determined by the matter density 
profile and the mixing angles. By using the Landau's method, the conversion probability Pjj for neutrinos to jump 
from one mass eigenstate to another at the high resonance layer can be expressed as [19I ] 



The factor F is given by 



Ph 



exp [— 7T7-F/2] — exp[— 7T7F/2sin 2 
1 - exp [-7r7F/2sin 2 6»] 



(9) 



F = < 



(l-tan 2 (9)7(l+tan 2 60 2 
(1 - tan 2 (9) 



2 E 

m=0 



/l/n-l\ 


1 
2 


V 2m J 


. m + 1 . 



(n e cx r) 
(n e oc r _1 ) 
(n e oc e~ r ) 

(tan26») 2m (n e cx r"), 



(10) 
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where n e is the electron density, r is the distance to the center of SN, and the adiabaticity parameter 7 is defined as 







Am 2 



7 



E\d\mi e /dr 



■sin20tan20, 



with Am 2 being the mass square difference of two mass eigenstates. 

For the SNs, n ps — 3, the expression of F in the case of n e oc r n in Eq. (|T0|) is 



m=0 



/ l/»-l\ 


1 

2 


\ 2m J 


TO + 1 



(tan20) 2 



In Eq. dHJl 



1/n-l \ _ (1/n-l)! 

2to J ~ (1/n-l -2to)!(2to)! : 



and 



with 



to + 1 



r/2 

J m = / (sin< 
'0 



(-1)' 



\2m, 



(7 m Jm-\-l 

' 



(2m- 1)!! 7T 
(2m)!! 2 



Eq. (|12[) can be expressed as a hypergeometric function: 



„ „ / n - 1 2n - 1 , „ 



In the case of e [0,7r/8], using the Euler integral representation [37[, one has 

r( c ) 



2 Fi(a,6;c;z) = m , 

r(6)r(c-6) 7 

We make the Taylor expansion for F near the point — = 0, 



t"- 1 (l-t) c - 6 - 1 (l-te)-°d*. 



F = F(0) + F\0)\^-j +F"(0)( ij +... + F( m )(0)^j +■■■ 
The first two coefficients in Eq. (fT8"|) can be obtained straightforwardly, 



F(0) = 1 - tan 2 0, F (0) = (1 - tan 2 



In (1 - tan 2 0) + 1 - 



1 + tan 2 fl 2 

* ln(l + tan 2 1 

tan 2 v 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 
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Comparing with the numerical result of the right-hand side of Eq. (jT7j) in the case i — ► 0, we find that the first two 
terms in Eq. (|18j) give dominant contributions and other items are negligible, so F can be approximately written as 

F = (1 - tan 2 6) (l - ^{M 1 - tan2 <?) + 1 - K 1 + tan2 #)/tan 2 0]ln(l + tan 2 6)}\ . (20) 

This expression is identical to that in [19I ]. 

For the case of n = —3, the comparison between the numerical result of the right-hand side of Eq. (fT7j) and that 
given by Eq. ([20]) is shown in Fig. 1. It can be seen that for n = — 3, Eq. ([20]) is a very good approximation to 2-P1 
in Eq. ([TTjl. 



B. SN shock wave effects and Ph 

The SN shock wave effects play an important role in the SN neutrino oscillations. As was pointed out in Ref. 



and further studied in Rcfs. [3i 



41 



, after the core bounce, the shock wave propagates inside the SN during the period 



of neutrino emission. It modifies the density profile of the star and the change is characterized by a density jump as 
shown in Fig. 2 [h| j4ll]. In several seconds, the shock wave may reach the resonance region where the conversion 
of different flavor SN neutrinos maximize, thus affecting the transition probability Pjj in the high resonance region 
~"^| 24 1 . The density distribution of SN might be divided into two phases roughly by time. Before the shock wave 



effects begin, i.e. t < Is, the SN matter density is [381 ]: 

poM^io--^) gfan*. (21) 

In Eq. ([2Tj) n = —2.4. We have checked that in this case Eq. ([20]) is a very good approximation to Eq. (fTTj) . The 
difference between the numerical results from these two equations are negligible when #13 = 8.8°. 
For t > Is, the shock wave effects set in and the matter density is given by 

p(r,t) = po(r)-r " (22) 

[l (r>r s ), 

where r s is the position of the shock wave front, f(r,t) is defined as [2] 

f(r,t) =exp{ [0.28- 0.69 ln(r s /fcm)][arcsin(l - r/r,)] 1 ' 1 }, (23) 
and £ is a typical ratio of the potential across the shock wave front, 

£ = V+/V. - 10, (24) 
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which measures the SN matter potential V(r) drop from 

V+ = lim V(r), (25) 

r— }r s 

to 

V- = lim V{r). (26) 

The SN matter potential is related to the SN electron density n e {r) by 

V(r) = V2G F n e (r) = V2G F N A p(r)Y e , (27) 

where Gf is the Fermi constant, Na is the Avogadro's number and Y e is the electron fraction. In the numerical 
calculations, we assume Y e = 0.5. In Eq. (22) and Eq. (23), a slightly accelerating shock-wave-front position r s is 
assumed with the explicit time dependence 
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r s (t) = -4.6 x lCT + 1.13 x 10 4 • t + 1 x 10 2 • t\ (28) 

where r s is in units of km and t in units of s. As shown in Fig. 2, there is a jump in the density curves when t > Is. 
Suppose on density curves there are three points (r = J^i,?^,^) satisfying the resonance condition, that is 

Am^- cos29 13 = 2V2G F n e (ri)E (i = 1,2,3), (29) 



where IAto^I |Ato3 2 | 2.4 x 10 _3 eF 2 [42jj. In this case, two of the points lie below the shock front and the third 
one beyond the shock front, i.e., n < r% < r s < r^, four crossing probabilities corresponding to these points are 
denoted by Phi, Phi, Ps and Ph3, respectively. We also define four densities: p res is the resonance density, which 
satisfies Eq. (|29|) [l|; p + is the density at r = r s corresponding to the matter potential V+ and p_ is the density at 
r = r s corresponding to the matter potential V- ; pb is the density at the bottom of the camber which is the minimum 
of the density profile below the shock front. Then in the case where three points satisfy the resonance condition, 
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Pb < Pres < P- , and Ph has the most complicated expression: 

Ph = (Phi + Pm + P s + Pm) - 2{P H iPm + PmPm + PhiPs + PmPm + PnnP. + PmP. 
+i{P H iP H 2Pm + PmP H 2Ps + PhiPhsPs + PbuPhsP.) - SPniPmPmPs 



+2(1- 


2P S - 


2P H3 + 


4,Ph 3 Ps)VPhiPh2(1 - 


Pm)(l- 


P H2 ) 


COS ( 


»12 


+2(1- 


2Phi 


-2P H3 


+ AP H1 P m )^P H2 P s {\ 


-Pm)(l 


-Ps) 


COS ( 


fo s 


+2(1- 


2P H 2 


-2P H3 


+ 4P H 2Ph 3 )\/PhiPs(1 


-Phi)(1 


-Ps) 


COS ( 


ha 


+2(1 - 


2P H i 


-2P H 2 


+ 4P ffl P ff2 ) v /P H3 P s (l 


-Pm)(l 


-Ps) 


COS ( 


h-3 


+2(1- 


2P S - 


2Phi + 


4PhiP s )VPh2Ph3{1- 


Pif 2 )(l - 


Pm) 


COS ( 


1>2D 


+2(1- 


2P S - 


2P H 2 + 


4Ph2P s WPhiPh3(1- 


Pm)(l- 


Pm) 


COS ( 


h 3 



(30) 



- 6 ^PmPm(l - Pm){l - Pm) cos <j> l2 ^P H3 P S {1 - Pm)(l - P s ) cos <p s3 , 

where the crossing probability Pm at r, (i = 1, 2, 3) can be calculated from Eq. (9), <f>ij = 1, 2, 3, s) is defined as 

« ^ ' dx^ ^ [Amii cos20 13 - 2£^(r)] 2 + (Am^ s in20 13 ) 2 - (31) 

The number of the points where the resonance condition occurs might be two, one, or even none. In the case where 
Pb < Pres < p+ and p res > p~, oi p res = Pb and p re s < P— , the resonance condition occurs at two points, r = ri,r% 
or r = ri, r 3 , then 

Pff = (Pm+Pffi+-P s )-2(PffiPffi+PffiP s +PffiP 8 )+4P ffl P ffi P s 



-2(1 - 2P s )y/P H1 P m (l - Pffi)(l - P ffl ) cos0 i; (32) 



+2(1 - 2P m ) VPfflP>(l - Pffi)(l - P-) cos^ is 

+2(1 - 2P H iWPhiPs(1 - Pm)(l - Ps) cos0 ls (1 = 2,3). 

When p res > or /9 res = /9& and p res > /?_, or p res < pi, and /? res < p_, the resonance condition occurs at only 
one point, r = r\ or r = r 3 , one obtains 

Ph = Puk + P s - 2P Hk P s + 2y/P Hk Ps{l - Pfffe)(l - P) cos</) fcs , (A; = 1,3). (33) 

If p_ < p res < pb, the resonance condition does not occur, it can be easily realized that 

Ph = Ps- (34) 

A summary for neutrino flavor conversions due to the neutrino shock wave effects and the MSW effects in various 
density regions is given in Table II [nil ]. 

Because the Daya Bay experimental result, #i 3 = 8.8° ± 0.8°, is quite different from the condition #i 3 < 3°, it is 
necessary to make clear the behavior of Ph for large #i 3 . From Eqs. (|9"))- (fTl"j) , it can be found that Ph depends on F 
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TABLE II: Summary of neutrino flavor conversions due to the neutrino shock wave effects and the MSW effects in various 
density regions. 



t < Is 


t > Is 




Pb > p- 


P 


- > Pb 




resonance in 


flavor conversion 


Resonance in 


flavor conversion 




region 


involved 


region 


involved 


Pm 


Pres > P+ 


Phi, Ps 


pres > P+ 


Phi, Ps 




P+ > Pres > Pb 


Pm, Pm, P s 


P+ > Pres > P- 


Pm, Pm, P s 




Pres — Pb 


Phi, Ps 


P- > Pres > Pb 


Pm, Pm, Ps, Pm 




Pb > Pres > P- 


Ps 


Pres Pb 


Pm, P s Pm 




P- > Pres 


Ps, Pm 


Pb > Pres 


Ps, Pm 



and 7 which is related to #13 and the neutrino energy E. 

Fig. 3, Fig. 4, and Fig. 5 illustrate the crossing probability Pjj as a function of the neutrino energy E, the time 
t and the mixing angle respectively. In Fig. 3, it can be seen that the value of Ph depends on the energy of 
SN neutrino. Whatever value #13 takes, the value of Ph has a great "jump" approximately at E = lOMeV. For 
013 = 3°, the curve of Pjj still has obvious continuous fluctuations from about 15MeV to higher energy. For 13 = 6° 
and #13 = 9°, the value of Ph changes smoothly and decreases slowly when E > 30MeV. Fig. 4 shows the curves of 
Ph for three typical neutrino energies when the time ranges from Os to 10s. We can see that as the energy increases 
the curve becomes fatter. In other words, the greater the neutrino energy, the longer time Ph keeps at high values. 
In general, the value of Ph reaches the maximum value when the time is between 4 — 6s. From Fig. 5(a), it is found 
that for a certain neutrino energy, at different times, the value of Ph changes smoothly in the range of #13 = 5° — 10°. 
However, the curve for t = 6s has rapid fluctuations between 0° and 5°. In Fig. 5(b), all the curves corresponding to 
three neutrino energies have obvious fluctuations when #13 is between 0° and 5°. This is far from the real #13 value. 
It can be seen from Figs. 3, 4, and 5 that Pjj is zero near the real value of 6*13(8.8°) when there are no shock wave 
effects. However, when the shock wave effects turn on Ph is not zero in a range of time. 



C. Collective effects and survival probability P vv 



The neutrino collective effects, which mechanism is totally different from the MSW effects, are caused by the 
neutrino-neutrino interaction inside the core-collapse SN. Recently, it has been realized to be a crucial feature at very 
high densities of the core and can change the emitted fluxes of different flavor neutrinos substanti ally 1221 12811 . 

Up to now, there have been a significant amount of studies on the neutrino collective effects 22] 28] 3JJ] 43j- 47 1. In 
Ref. [3], it was shown that the collective effects depend on the inherent features of SN neutrinos, such as their initial 
total energy, relative luminosities of different flavors, and the neutrino mass hierarchy [l| {3^] . 

In order to study the collective effects quantitatively, we set P vu as the survival probability that the neutrinos 
(antincutrinos) v(p) remain as their original states after the collective effects. In Ref. the authors introduced a 
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simplified picture to describe the characteristics of the collective effects 

1 [E<E C ), 



(35) 
(E>E C ), 



for neutrinos and 



(36) 

for antineutrinos, where Ec is a critical energy. We take Ec = lOMeV in our later calculation (22 



D. Earth matter effects 

In this subsection, we briefly review known analytical results about the Earth matter effects, which were studied 
in Rcf. [sj] for the sake of completeness and self-consistency of this paper. If a neutrino reaches the detector with an 
incident angle a (see Fig. 1 in Ref. [lo|]), the distance that the neutrino travels through the Earth to the detector, L, 
and the distance of the neutrino to the center of the earth, x, can be given by 



L = (-R E + h) cos a + \J R E - (Re - h) 2 shr a, 

x = \/(-Re + h) 2 + (L — x) 2 + 2(R E - h)(L - x) cos a, 

where h is the depth of the detector in the Earth, x the distance that the neutrino travels into the Earth, and Re the 
radius of the Earth. 

Let Pi e be the probability that a neutrino mass eigenstate Vi enters the surface of the Earth and arrives at the 
detector as an electron neutrino v e , then we have [49 1 



P 2e = sin 2 6ia + \ sin 2 20 12 f ' dxV E (x) sin , (37) 

1 Jxn 



where Ve(x) is the potential v e experiences due to the matter density pe(x) inside the Earth 



V E (x) = V2G F N APE (x)Y e , (38) 



and f/'^;, is defined as 



b 

dxA m (x), 



a 



where 

A m (x) = ^k,J (cos26 12 - e(x)) 2 + sin 2 20 12 , (39) 
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with #12 12 34.5°, Ato|x = 7.7 x 10~ 5 eT^ 2 , e(x) = 2£ , V E (x)/Am| 1 |42j. We will use the realistic matter density profile 
inside the Earth to compute the Earth matter effects. 

IV. DETECTING SN NEUTRINOS ON THE EARTH 

As mentioned in Rcfs. fllltj, there are three reaction channels with which one may detect SN neutrinos at the 
Daya Bay experiment. In this section, we review the three reactions and the calculation methods of event numbers. 
The reaction formulas and the means of calculation are still applicable when #13 is large. 

A. Three reactions 

A detailed description of the Daya Bay experiment can be found in Rcfs. There are eight detectors located 

at different sites of the Daya Bay experiment. The total detector mass is about 300 tons and the depth of the detectors 
h i2 400 m. The Daya Bay Collaboration uses Linear Alkyl Benzene (LAB) as the main part of the liquid scintillator. 
LAB has a chemical composition including C and H with the ratio of the number of C and H being about 0.6. Then, 
the total numbers of target protons, electrons, and 12 C are 

= 2.20 x 10 31 , N^ e) = l.oi x 10 32 , = 1.32 x 10 31 . 

In the Daya Bay experiment, there are three reactions which can be used to detect SN neutrinos: the inverse beta 
decay^neutrino-electron reactions, and neutrino-carbon reactions. Their effective cross sections are given as follows 

(1) The cross section for the inverse beta-decay is 

a{D eP ) = 9.5 x \Q- AA {E{MeV) - 1.29) 2 cm 2 , (40) 

where the reaction threshold is E t h = 1.80MeV. 

(2) The neutrino-electron scattering interactions can be divided into changed current and neutral current ones. The 
total cross sections have the following forms: 

(cr(v e e -> v e e)) = 9.2 x 10~ 45 £(M eV) 2 cm 2 , 
(a(D e e -> v e e)) = 3.83 x 10~ 45 E(MeV) 2 cm 2 , 
(cr(iV, T e -> v^ T e)) = 1.57 x 10" 45 ^(Af eV) 2 cm 2 , 

(a(u M , T e -> 2^, T e)) = 1.29 x 10~ 45 ^(Afel/) 2 cm 2 . (41) 

(3) The cross sections for the neutrino-carbon reactions are 

(<j( 12 C(i/ e ,e-) 12 A0) = 1.85 x 10- 43 £(AfeT0 2 cTO 2 , 

(aC 2 C(9 e ,e+) 12 B)) = 1.87 x IQ- 42 E(MeV) 2 cm 2 , (42) 
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for the charged-current capture, and 

{(j{vl 2 C)) = 1.33 x W- i3 E(MeV) 2 cm 2 , 

{a{vl 2 C)) = 6.88 x lO- 43 E(MeV) 2 cm 2 , (43) 
(<?(v^t(v^t) 12 C)) = 3.73 x 10~ 42 E(MeV) 2 cm 2 , 

for the neutral-current capture. 

We should note that when neutrino oscillations are taken into account, the oscillations of higher energy v x into 
v e result in an increased event rate since the expected v e energies in the absence of oscillations are just at or 
below the charged-current reaction threshold. This leads to an increase by a factor of 35 for the cross section 
(cr( 12 C(V e , e~) l2 N)) if we average it over a v e distribution with T = 8MeV rather than 3.5MeV. Similarly, the cross 
section (er( 12 C(j> e , e + ) 12 B)) is increased by a factor of 5. 

In the following subsections, we will discuss the SN neutrino fluxes taking into account the four physical effects and 
calculate event numbers. 

B. Calculation of event numbers 

With all of the four physical effects being taken into account, the SN neutrino fluxes at the detector are expressed 

as 

F? = pF^ + (l-p)F^, 
2F% = (l-p)F^ + (l+p)F^, 

2F° = (1 - p)F^ + (1 + p)F^ , (44) 
where the survival probabilities p and p are given by 

P = P2e[P H Puu + (l-P H )(l-P,u)}, 

p = (l-P 2e )P uu , (45) 

for the normal mass hierarchy and 

P = PlePvv, 

p = (1 - P2e)[P H Puu + (1 - P H )(1 - Puu)}, (46) 

for the inverted mass hierarchy. 

Therefore, the event numbers N(i) of SN neutrinos in the reaction channel 'V can be calculated by integrating 
over the neutrino energy, the product of the target number Nt, the cross section of the given channel cr(i), and the 
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neutrino flux function at the detector F® , 

N(i)=N T j &E-a(i)-^-FS, (47) 
where a stands for the neutrino or antineutrino of a given flavor, and D is the distance between the SN and the Earth. 

C. Scanning over the relevant parameters 

In this subsection, we scan the relevant parameters in the two parametrization forms for neutrino energy distribution. 
In Section II, we gave the ranges of the neutrino temperatures and pinching parameters of the spectra in the Fcrmi- 
Dirac distribution, as listed in Eq. (j4|), and the ranges of the average energies of neutrinos and dimensionlcss pinching 
parameters in the parametrization of SN neutrino fluxes given by the Garching group, as listed in Eq. (|7J|. It is 
expected to obtain the maximum and minimum values of neutrino event numbers in the Daya Bay experiment from 
our calculation results. To achieve this objective, scanning over the ranges of all parameters related to the calculation 
of four physical effects on detecting SN neutrinos is necessary. Notice that the luminosity ratios of different flavor 
neutrinos in Eq. ([5]) and Eq. ([SJ should be considered as well. In fact, as will be shown in the next subsection, the 
luminosity ratios do have effects on neutrino event numbers. 

Comparing the three reaction channels in Eqs. (|40 |) -(|43 |) . it can be seen that the cross sections for the neutrino- 
electron scattering channel are much smaller than the other two reaction channels Hence, we will only consider 
the inverse beta-decay and the neutrino-carbon reactions in the following analysis. From Eq. (|40l) . it can be seen that 
the inverse beta-decay does not involve any parameters about v e since v e is not involved in the inverse beta-decay. 
Based on simulation results the luminosity of v e and v e can be taken to be equal , so we can define 



ha. = ±2± = J_ (48) 

Scanning over the ranges of all the parameters, we have numerical results for the event numbers. For the Fermi-Dirac 
distribution in Eq. ((3]), we find the following two groups of parameters correspond to the maximum and minimum 
event numbers, respectively, 

{Max) T Pc = 5MeV, T„ x = 9MeV, rj Pe = 2, r} Vx = 2, M = 2, (49) 



(Min) T Pe = 5MeV, T Vx = 7MeV, t]^ = 2, ?7„ x = 0, M = 0.5, (50) 
for the inverse beta-decay, and 



(Max) T Uc = 3MeV, T Pc = 5MeV, T Vx = 7MeV, r] Ve = 3, 77^ = 2, r] Vx = 0, M = 2, 



(51) 
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(Min) T Uc = AMeV, T De = 6MeV, T Vx = 9MeV, n Ve = 5, n„ e = 2.5, n Vx = 2, M = 0.5, (52) 

for the neutrino-carbon reactions. 

From Eq. ([5]) it can be seen that M varies between two extreme values, 0.5 and 2. In order to see the influence of 
the luminosity ratio itself on event numbers, we only change the values of M in Eqs. (|4U]) - (|52"|) to the other extreme 
values, respectively, while keeping other parameters unchanged, then we obtain the following comparison groups of 
parameters: 

(Max - C) T Pc = 5MeV, T„ x = 9MeV, n Pe = 2, r\ Vx =2, M = 0.5, (53) 

(Min - C) T Pe = 5MeV, T Vx = 7MeV, r toe = 2, n Vx =0, M = 2, (54) 
for the inverse beta-decay reactions, and 

(Max - C) T Ve = 3MeV, = 5MeV, T Vx = 7MeV, r) Ve = 3, n p<! = 2, n Vx =0, M = 0.5, (55) 

(Min - C) T Ve = iMeV, = 6MeV, T Vx = 9MeV, n Ve = 5, n 9c = 2.5, r\ Vx = 2, M = 2, (56) 

for the neutrino-carbon reactions. 

In the same way, after scanning over the ranges of all the parameters in Eqs. ([7]) and ([8|) we obtain the following 
parameter values for the distribution in Eq. ((6]): 

(Max) (E Pe ) = UMeV, (E Vm ) = l8MeV, /3 Pc = 3.5, p Va = 3.5, M = 2, (57) 

(Min) (Ec c ) = UMeV, (E Vx ) = !5MeV, /3 Pe = 3.5, /3„ x = 6, M = 1.25, (58) 
for the inverse beta decay, and 

(Max) (E Ue ) = UMeV, (E P J = UMeV, (E V J) = VbMeV, /3 Ue = 3.5, /3 Pa = 3.5, /3 Vx = 6, M = 2, (59) 

(Mm) (E Ve ) = lbMeV, (Ep.) = 15MeV, (E Vm ) = 18MeV, p Ue = 3.5, /3 Pc = 3.5, p Vm = 3.5, M = 1.25, (60) 

for the neutrino-carbon reactions. 

Similar to Eqs. (|53|) - (|56[) . we have the comparison groups of parameters for the distribution in Eq. ([5]): 



(Max - C) (Eb.) = UMeV, (E„ m ) = 18MeV, fr c = 3.5, = 3.5, M = 1.25, 



(61) 
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(Min - C) (E Ve ) = 12MeV, (E v .) = 15MeV, [3^ = 3.5, =6, M — 2, (62) 
for the inverse beta decay, and 

(Max - C) (E Ue ) = 12MeV, (E s .) = 12MeV, (E Vx ) = 15MeV, p Ue = 3.5, fo B = 3.5, =6, M = 1.25, (63) 

(Min - C) (E Vb ) = 15MeV, (E Pe ) = 15MeV, (E V J = l&MeV, = 3.5, /3 Pc = 3.5, P„ m = 3.5, M = 2, (64) 

for the neutrino-carbon reactions. 

It is noted that no matter how the value of M changes, the total energy of all flavor neutrinos is a constant and 
the results with parameters in comparison groups are always between the maximum and minimum event numbers. 

D. The SN neutrino event numbers under the influence of four physical effects 

In this subsection, we give the numerical results of SN neutrino event numbers detected at the Daya Bay experiment. 
Consider a "standard" supernova at a distance D = lOkpc from the Earth, which releases total energy Eb = 



3 x 10 58 erg, and take r = 3s as the decay time of its luminosity [30|-[32j. The values of relevant parameters have 
already been given in detail in previous sections. Given the Daya Bay experimental result, we take $13 = 8.8° in our 
calculations. 

Firstly, we calculate the neutrino event numbers with the neutrino spectra of the Fermi-Dirac distribution. With 
the influence of all the four physical effects being taken into account, the neutrino event numbers TV are plotted as 
a function of the incident angle of the neutrino a in Fig. 6 and Fig. 7 for the inverse beta-decay and the neutrino- 
carbon interactions, respectively. In order to show the influence of the luminosity ratio, we show the results with 
"Max (Min)" group of parameters and those with "Max — C (Min — C)" group of parameters in same figures. For 
example, in Fig. 6(a) "Max" represents the results with parameters in Eq. (|49|) and "Max — C" represents those 
with parameters in Eq. (|53|) . It is found that event numbers depend on the mass hierarchy and parameter values. 
The maximum variation of neutrino event numbers appears at a ~ 93° when a changes due to the Earth matter 
effects. The variations in the inverse beta-decay are more obvious than those in the neutrino-carbon reactions. 

In Table. Ill and Table. IV, we sum up the neutrino event numbers for the two reaction channels with the Fermi- 
Dirac parametrization. The numerical results show that the event numbers and the change rates due to the Earth 
matter effects depend on the parameters T a and rj a . This was observed previously in Ref. [9(. It is found in the 
present work that when the values of T a and rj a keep unchanged, the luminosity ratio M plays an important role 
in determining the event numbers. For instance, comparing the results with the "Max" and "Max — C" groups of 
parameters for the inverse beta-decay, in the case of normal mass hierarchy, when M = 0.5 the event numbers at the 
incipient angle and the change rate are 110.21 and 4.81%, respectively; while when M = 2 these two numbers are 
132.0 and 14.08%, respectively, which are much larger. 

Next, we discuss the results with the "Garching" distribution given in Eq. ([6]) with the parameter values being 
listed in Eqs. (|57]l- ([64]) . Similar to Fig. 6 and Fig. 7, Fig. 8 and Fig. 9 show the behavior of event numbers for 



the inverse beta-decay and the neutrino-carbon interactions, respectively. The summary of event numbers for the two 
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reaction channels with the "Garching" distribution is given in Table. V and Table. VI. From these figures and tables, 
we also find that the event numbers and the change rates due to the Earth matter effects depend on the parameters 
(E a ) and (3 a . Furthermore, the luminosity ratio M has an important influence on the event numbers and the change 
rates while keeping other parameters unchanged. For example, comparing the results with the "Max" and "Max — C" 
groups of parameters for the neutrino-carbon interactions, in the case of normal mass hierarchy, when M = 1.25 the 
event numbers at the incipient angle and the change rates are 61.49 and 0.15%, respectively; while when M = 2 these 
two numbers are 77.4 and 0.22%, respectively, which are much larger. 

Comparing the results with two parametrization forms of neutrino energy distribution, we can easily see that for 
the inverse beta-decay, event numbers with the Fermi-Dirac distribution are larger than those with the "Garching" 
distribution. However, for the neutrino-carbon interactions, the event numbers with the Fermi-Dirac distribution are 
smaller than those with the "Garching" distribution. On the other hand, the change rates of event numbers in the 
two reaction channels due to the Earth matter effects in the case of the Fermi-Dirac distribution are always larger 
than those in the case of the " Garching" distribution. 

Based on the results in Tables III- VI, in Table VII we give neutrino event numbers detected at the Daya Bay 
experiment when all the uncertainties are taken into account in the cases of normal and inverted mass hierarchies, 
respectively. We can see that the event numbers range from 63.66 ~ 243.51 and 16.12 ~ 94.9 for the inverse beta-decay 
and the neutrino-carbon interactions, respectively. 



V. SUMMERY AND DISCUSSIONS 



Given the new experimental result about #13 from the Daya Bay Collaboration, we deduce the expression of the 
neutrino conversion probability in the high resonance region inside SN, Pjj, in the case of large #13 (i2 8.8°), by 
applying the available knowledge for Ph- Ph is expressed in the form of hypcrgcomctric function. In the derivation, 
we take the shock wave effects into account. We give numerical results of Ph as functions of #13, t, and E. It is found 
that Ph is zero near the real value of #13 when there are no shock wave effects. However, it is not zero in a certain 
region of time (roughly 3s ~ 8s depending on neutino energies) when shock wave effects are taken into account. Our 
work is different from previous studies which were usually done in the case of small #13 (< 3°) f^fiol]. 

We consider the influence of all the four physical effects on the detection of SN neutrinos, including the MSW 
effects, the SN shock wave effects, the neutrino collective effects, and the Earth matter effects. Scanning over all the 
relevant parameters in the two different parametrization forms of neutrino energy distribution, we calculate the event 
numbers for two reaction channels, the inverse beta-decay and the neutrino-carbon reactions, both of which can be 
detected at the Daya Bay experiment. It is found that the event numbers depend on the parameters T, A , //,,. (E a ), 
/3 a , and L a , as well as the mass hierarchy. Wc emphasize that the event numbers depend on the luminosity ratio 
substantially. We give the range of SN neutrino event numbers detected at the Daya Bay experiment. 

Although a lot of effort has been made on identifying the four physical effects on detection of SN neutrinos, there 
are still a lot of problems which need to be solved. One example is the neutrino collective effects in the case of large 
#13. Progress in this direction will affect the detection of SN neutrinos. Now the Daya Bay II experiment is under 
consideration. Its detector mass will be about 70 times of the total detector mass of the Daya Bay experiment. This 
will make it much more possible to detect SN neutrinos in the future. 
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Figure Captions 

Fig. 1 Numerical result of F as a function of 9 for n = — 3. The solid and dashed curves represent the results when 
F takes the expression of Eq. (fl7j) and Eq. ([20]) , respectively. 



Fig. 3 The crossing probability at the high resonance region Ph as a function of the neutrino energy E for three 
neutrino mixing angles at t = 6s. The solid, dashed and dotted curves correspond to #13 = 3°, 6°, 9°, respectively. 

Fig. 4 The crossing probability at the high resonance region Ph as a function of time t for three neutrino energies at 
#13 = 9°. The solid, dashed and dotted curves correspond to neutrino energy E = 11, 16, 25MeV, respectively. 

Fig. 5 The crossing probability at the high resonance region Pjj as a function of the neutrino mixing angle $13: 
(a) for three different times at E = HMeV. The solid, dashed and dotted curves correspond to t = 2s, 4s, 6s, 
respectively; (b) for three different neutrino energies at t = 6s, The solid, dashed, and dotted curves correspond to 
E = 11, 16, 25MeV, respectively 

Fig. 6 The event numbers for the inverse beta-decay interaction with the parametrization form in Eq. (J3]). "a" is the 
incident angle of the SN neutrino reaching the detector. "N(I)" represents normal (inverted) mass hierarchy. In Fig. 
6(a) "Max" and "Max-C" correspond to parameter values listed in Eq. ([49]) and Eq. ([53]) . respectively; in Fig. 6(b) 
"Min" and "Min-C" correspond to parameter values listed in Eq. (|50|) and Eq. ([54]) . respectively. 

Fig. 7 The event numbers for the neutrino-carbon interactions with the parametrization form in Eq. (|3|). "a" is the 
incident angle of the SN neutrino reaching the detector. "N(I)" represents normal (inverted) mass hierarchy. In Fig. 
7(a) "Max" and "Max-C" correspond to parameter values listed in Eq. ([51]) and Eq. (f55|) . respectively; in Fig. 7(b) 
"Min" and "Min-C" correspond to parameter values listed in Eq. ([52]) and Eq. (|56|) . respectively. 

Fig. 8 The event numbers for the inverse beta-decay interaction with the parametrization form in Eq. (j6]). "a" is the 
incident angle of the SN neutrino reaching the detector. "N(I)" represents normal (inverted) mass hierarchy In Fig. 
8(a) "Max" and "Max-C" correspond to parameter values listed in Eq. ([57)) and Eq. ([6Tj) , respectively; in Fig. 8(b) 
"Min" and "Min-C" correspond to parameter values listed in Eq. ([58]) and Eq. ([62]) . respectively. 

Fig. 9 The event numbers for the neutrino-carbon interactions with the parametrization form in Eq. (j6]). "a" is the 
incident angle of the SN neutrino reaching the detector. "N(I)" represents normal (inverted) mass hierarchy. In Fig. 
9(a) "Max" and "Max-C" correspond to parameter values listed in Eq. ([59]) and Eq. ([63]) . respectively; in Fig. 9(b) 
"Min" and "Min-C" correspond to parameter values listed in Eq. ([60]) and Eq. ([64]) . respectively. 



Fig. 2 The changes of the density of the supernova under the influence of the shock wave effects 
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TABLE III: Summary of event numbers of the inverse beta-decay interaction for the parametrization form given in Eq. ([3|. 
"Max" and "Max — C" correspond to parameter values in Eq. (|49[) and Eq. I|53[l. respectively; "Min" and " Min — C" 
correspond to parameter values in Eq. (|50[) and Eq. (|54[) , respectively. "N(I)" represents normal (inverted) mass hierarchy. 
The numbers in the columns "Incipient" and "Min" are the event numbers when the SN neutrino incident angle is zero and 
is the angle in the column "Angle", respectively. The column "Angle" gives the angles at which the event numbers are the 
minimum and the Earth matter effects are the strongest. "Ratio" gives the percentages of the Earth matter effects. 



Reaction 


Conditions 


Hierarchy 


Incipient 


Min 


Angle 


Ratio 




Max 


N 


132.00 


113.41 


94 


14.08% 






I 


245.51 


243.13 


94 


0.97% 




Max - C 


N 


110.21 


104.91 


94 


4.81% 


Ve +P 




I 


107.84 


107.22 


94 


0.57% 




Min 


N 


100.33 


99.67 


94 


0.66% 






I 


79.26 


79.21 


94 


0.06% 




Min-C 


N 


107.29 


99.37 


94 


7.38% 






I 


174.07 


172.94 


94 


0.65% 



TABLE IV: Summary of event numbers of the neutrino-carbon interactions for the parametrization form given in Eq. ([3|. 
"Max" and "Max — C" correspond to parameter values in Eq. (|51l) and Eq. (|55[) . respectively; "Min" and "Min — C" 
correspond to parameter values in Eq. I|52|) and Eq. (|56p . respectively. "N(I)" represents normal (inverted) mass hierarchy. 
The numbers in the columns "Incipient" and "Min" are the event numbers when the SN neutrino incident angle is zero and 
is the angle in the column "Angle", respectively. The column "Angle" gives the angles at which the event numbers are the 
minimum and the Earth matter effects are the strongest. "Ratio" gives the percentages of the Earth matter effects. 



Reaction 


Conditions 


Hierarchy 


Incipient 


Min 


Angle 


Ratio 




Max 


N 


50.86 


50.68 


92 


0.35% 






I 


63.01 


62.92 


92 


0.14% 




Max - C 


N 


23.79 


23.76 


93 


0.13% 






I 


26.97 


26.95 


92 


0.07% 




Min 


N 


16.15 


16.11 


93 


0.25% 






I 


18.19 


18.17 


93 


0.11% 




Min - C 


N 


33.84 


33.62 


93 


0.65% 






I 


42.27 


42.18 


93 


0.21% 
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TABLE V: Summary of event numbers of the inverse beta-decay interaction for the parametrization form given in Eq. (JSJ - 
"Max" and "Max — C" correspond to parameter values in Eq. (|57[) and Eq. I|61[l. respectively; "Min" and " Min — C" 
correspond to parameter values in Eq. (|58[) an d Eq. (|62[) , respectively. "N(I)" represents normal (inverted) mass hierarchy. 
The numbers in the columns "Incipient" and "Min" are the event numbers when the SN neutrino incident angle is zero and 
is the angle in the column "Angle", respectively. The column "Angle" gives the angles at which the event numbers are the 
minimum and the Earth matter effecs are the strongest. "Ratio" gives the percentages of the Earth matter effects. 



Reaction 


Conditions 


Hierarchy 


Incipient 


Min 


Angle 


Ratio 




Max 


N 


78.18 


73.21 


93 


6.36% 






I 


135.58 


134.81 


92 


0.57% 




Max - C 


N 


74.59 


71.13 


93 


4.64% 


i> e +p 




I 


106.89 


106.38 


93 


0.48% 




Min 


N 


64.77 


63.92 


92 


1.31% 






I 


78.98 


78.84 


92 


0.18% 




Min - C 


N 


65.41 


63.66 


92 


2.68% 






I 


99.31 


99.01 


92 


0.30% 



TABLE VI: Summary of event numbers of the neutrino-carbon interaction for the parametrization form given in Eq. ||6j. 
"Max" and "Max — C" correspond to parameter values in Eq. 1)591) and Eq. I)63p . respectively; "Min" and "Min — C" 
correspond to parameter values in Eq. (|60[) and Eq. (|64[) . respectively. "N(I)" represents normal (inverted) mass hierarchy. 
The numbers in the columns "Incipient" and "Min" are the event numbers when the SN neutrino incident angle is zero and 
is the angle in the column "Angle", respectively. The column "Angle" gives the angles at which the event numbers are the 
minimum and the Earth matter effects are the strongest. "Ratio" gives the percentages of the Earth matter effects. 



Reaction 


Conditions 


Hierarchy 


Incipient 


Min 


Angle 


Ratio 




Max 


N 


77.40 


77.23 


95 


0.22% 






I 


94.90 


94.84 


92 


0.06% 




MI ax - C 


N 


61.49 


61.40 


92 


0.15% 






I 


73.94 


73.90 


92 


0.05% 




Min 


N 


50.16 


50.05 


92 


0.22% 






I 


60.52 


60.47 


92 


0.08% 




Min - C 


N 


63.20 


63.03 


92 


0.27% 






I 


77.70 


77.62 


92 


0.10% 



TABLE VII: The event number ranges in the Daya Bay experiment with all the uncertainties taken into account. "N(I)" 
represents normal (inverted) mass hierarchy. 



Reaction 


Hierarchy 


Max 


Min 


Range 


V e +P 


N 


132.00 


63.66 


63 - 246 


I 


245.51 


78.84 


v+ 12 C 


N 


63.20 


16.11 


16 ~ 95 


I 


94.90 


18.17 



